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ABSTRACT 

^ ! We are conducting a survey of Galactic shell-type supernova remnants (SNRs) known or suspected to possess non- 
| thermal components to their X-ray emission using new and archived observations made with such X-ray satellites as 
"^t" ' ROSAT, ASCA, RXTE, Chandra and XMM — Newton. This research is intended to probe the phenomenon of 
q . cosmic-ray acceleration by Galactic SNRs and estimate the maximum energy of cosmic-ray electrons accelerated by 
^ | these sources. To illustrate this work, we examine the X-ray spectrum of the northwestern rim of an SNR suspected 
O ■ to have a non-thermal component to its X-ray emission, G156.2+5.7 (RX J04591+5147), over the energy range of 
^ . «0.7-12.0 keV using observations made by the ASCA GIS and the RXTE PC A. We compare fits made to the non- 
thermal component using two models, a simple power law and SRCUT. Both models give acceptable fits: the photon 



O ' index derived from the fit made with the power law model (T = 2.0_q' 5 ) is comparable to values obtained for the 
£j bright rims of other SNRs with hard X-ray spectra. Using the SRCUT model, we derive a value of 2.42^23 X 10 17 
j3 Hz for the cutoff frequency z^cutoff : based on this value and assuming a mean magnetic field strength of 14 /j,G, we 
estimate the cutoff energy E cnto g of cosmic -ray electrons accelerated by G 156. 2+5. 7 to be rs 32 TeV. This energy 
'"j ■ value is well short of the "knee" feature of the cosmic -ray spectrum. 



INTRODUCTION 



In recent years the study of non-thermal X-ray emission from Galactic supernova remnants (SNRs) has attracted 
a considerable amount of interest within the astronomical community. Extensive research in this field commenced 
with the discovery that the X-ray emission from the Galactic SNR SN 1006 is predominately non-thermal (Koyama 
et al. 1995, Allen et al. 2001). Since this discovery, the number of Galactic SNRs which are known to possess spectra 
dominated by non-thermal X-ray emission has steadily increased and currently includes G347. 3-0.5 (Koyama et al. 
1997, Slane et al. 1999, Pannuti et al. 2003), G266.2-1.2 (Slane et al. 2001) and G28.6-0.1 (Bamba et al. 2001, 
Koyama et al. 2001) as well as SN 1006. Generally speaking, all of these Galactic SNRs are X-ray luminous, radio 
weak sources with shell-type morphologies expanding into low density ambient media. It is generally believed that 
the non-thermal X-ray emission is in fact synchrotron emission from electrons accelerated along the shock front of 
the SNR (e.g., Reynolds 1996, 1998; Reynolds and Keohane 1999), although others have argued that for at least some 
SNRs, the emission is instead produced by non-thermal bremsstrahlung (e.g. Vink and Laming 2002). 

If the observed non-thermal emission does in fact have a synchrotron origin, it is possible to estimate the max- 
imum energy of cosmic-ray electrons accelerated by these SNRs through an analysis of this emission. While it 
is currently believed that Galactic SNRs accelerate cosmic-ray particles to the well-known "knee" feature of the 
cosmic-ray spectrum (i^knee ~ 3000 TeV), only a scant amount of data is available to test this hypothesis. To probe 
this issue in more detail, we are conducting an X-ray survey of a sample of Galactic SNRs using new and archived 



Table 1 . ASCA and RXTE Observations of the Northwestern Rim of G1 56.2+5.7 (RX J04591 +51 47) 



R.A. Decl. 

Satellite and SeqID or Observation (h m s) (° m s) 

Instrument ObsID Date (J2000.0) (J2000.0) 



Radius of 
Field of View 
(arcmin) 



Effective Range of 
Exposure Sampled 



Time Energies 
(Seconds) (keV) 



ASCA GIS2+GIS3 50011000 1993 Sep 25 04 59 14 +52 29 56 
ASCA GIS2+GIS3 57059000 1999 Mar 5 04 55 36 +52 06 52 
RXTEPCA 50150-01 2002 Jan 25 04 56 45 +52 21 00 



22 
22 
60 



32964 0.7-2.5 
42570 0.7-2.5 
16160 3.0-12.0 



observations made with such satellites as the Rontgensatellit (ROSAT), the Advanced Satellite for Cosmology and 
Astrophysics (ASCA), the Rossi X-ray Timing Explorer (RXTE), the X-ray Multiple Mirror Mission (XMM-Newton) 
and the Chandra X-ray Observatory. We are extracting emission spectra from the X-ray bright rims of Galactic SNRs 
known or suspected to possess non-thermal components to their X-ray emission. To illustrate this work, in this pa- 
per we consider one of the sources in our sample, G156.2+5.7 (RX J04591+5147): this SNR is suspected to possess 
both thermal and non-thermal components to its X-ray emission (Yamauchi et al. 1993, Yamauchi et al. 1999, Tomida 
1999). We analyze observations made of the northwestern rim of this SNR using the Gas Imaging Spectrometer (GIS) 
aboard ASCA and the Proportional Counter Array (PCA) aboard RXTE. A detailed summary of these observations 
is provided in Table 1. In Figure 1, we present an X-ray image of G156.2+5.7 that was prepared using six pointed 
observations made with the Position Sensitive Proportional Counter (PSPC) aboard ROSAT. 



THE GALACTIC SNR G156.2+5.7 (RX J04591+5147) 

The Galactic X-ray SNR G156.2+5.7 (RX J04591+5147) was discovered during the ROSAT All-Sky Survey 
and first described by Pfeffermann et al. (1991). Using a Sedov model to fit its X-ray spectrum, those authors derived 
an age of 26000 yr for G156.2+5.7, a distance of 3 kpc (which is supported by an independent distance estimate of 
1-3 kpc based on the column density toward this source as well as HI observations) and an extremely low ambient 
density of 0.01 cm -3 . Assuming a distance of 3 kpc, the X-ray flux from this SNR given by Pfeffermann et al. (1991) 
(1.9x 10 -10 ergs cm~ 2 sec -1 between 0.1 and 2.4 keV) indicates a luminosity of 2.0 x 10 35 ergs sec -1 , placing this 
source among the ten most luminous known Galactic X-ray SNRs. Complementary radio observations of G156.2+5.7 
presented by Reich et al. (1992) indicated a modest flux density of 4.2+0.1 Jy at 1 GHz and an extremely low surface 
brightness of 5.8 x 10~ 20 ergs cm~ 2 sr _1 . Reich et al. (1992) also noted that the radio emission from the rims of 
this SNR is highly polarized, consistent with a synchrotron origin. These characteristics - high X-ray luminosity 
but low radio luminosity, low ambient medium density, and polarized radio emission from the rims of the SNR - 
are generally in common with the other Galactic SNRs known to possess X-ray spectra dominated by non-thermal 
emission. Additional X-ray observations of G156.2+5.7 have been presented using observations made with ROSAT, 
Ginga (Yamauchi et al. 1993) and ASCA (Yamauchi et al. 1999). Both of those papers commented on the harder 
component to the X-ray spectrum of this SNR, with Yamauchi et al. (1999) deriving a photon index of r=2.llQ'| for 
a power-law fit to this component over the ASCA energy band. This index is comparable to the indicies derived by 
power-law fits to the hard components of the ASCA X-ray spectra from other Galactic SNRs. 

OBSERVATIONS OF G156.2+5.7 

ASCA GIS2+GIS3 Observations 

Both the GIS and the Solid-State Imaging Spectrometer (SIS) aboard ASCA (Tanaka et al. 1994) observed two 
different positions along the northwestern rim of G156.2+5.7, but becuase the GIS has a larger field of view than 
the SIS we only consider data from that instrument in this work. The GIS consists of two units (denoted as GIS2 
and GIS3) of imaging gas scintillation proportional counters with a sealed-off gas cell equipped with an imaging 
phototube. The GIS is sensitive over the energy range of approximately 0.7 - 10 keV and has an energy resolution of 
AE/E = 0.078 at 6 keV (Ohashi et al. 1996, Makishima et al. 1996). The reduction process filtered data for elevation 
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Fig. 1 . ROSAT PSPC X-ray image of G1 56.2+5.7 
depicting emission over the energy range of 0.4- 
2.0 keV (courtesy of Steven Snowden). The two 
small circles represent the fields of view of the 
ASCA GIS observations, while the large circle 
represents the field of view of the RXTE PCA 
observation (see Table 1 and "Observations of 
G 156.2+5.7.") 
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Fig. 2. X-ray spectra of the northwestern rim of 
G1 56.2+5.7 as observed by the ASCA GIS and 
the RXTE PCA (see Table 1 and "Results"). 



angle, stable pointing directions, the South Atlantic Anomaly and cut-off rigidity. The data were also screened for 
characteristics of the GIS internal background and screening based on event locations and rise time. The standard 
script "ASCASCREEN" was run to accomplish this reduction. To boost the signal-to-noise, for both pointings we 
combined the spectra observed by the GIS2 and GIS 3 units and considered these combined spectra in our analysis. 

RXTE PCA Observations 

Because of the rather large angular extent of G156.2+5.7 (sa 1 10 arcminutes, Green 2001), it was not possible to 
observe the entire SNR with a single pointing of the PCA, so we instead observed just the luminous northwestern rim. 
The PCA is a spectrophotometer comprised of an array of five coaligned proportional counter units that is sensitive 
to photons that have energies between approximately 2 and 60 keV. The energy resolution AE/E of the array is 0.18 
at 6 keV, and the maximum on-axis collecting area is about 6000 cm 2 at 9.72 keV. The PCA data were screened to 
remove the time intervals during which (1) one or more of the five proportional counter units is off, (2) G156.2+5.7 is 
less than 10° above the limb of the Earth, (3) the background model is not well-defined and (4) the pointing direction 
of the detectors is more that 0°.02 from the nominal pointing direction in either right ascension or declination. 

RESULTS 

We first considered the two ASCA GIS observations of the northwestern rim: these observations were sensitive 
to thermal emission from G156. 2+5.7 and we used the data from these observations to model this type of emission. 
We point out that the point source AX J0500+5238 lies along the northwestern rim and is believed to be a background 
galaxy unrelated to G156.2+5.7 (Tomida 1999, Yamauchi et al. 1999). We extracted spectra for two regions along 
this rim as sampled by the two ASCA GIS pointings (excluding flux from AX J0500+5238) and fit the spectra using 
a thermal RAYMOND model (Raymond and Smith 1977) with abundances frozen to solar values. Photoelectric 
absorption along the line of sight was fit using the WABS model, which is based on the Wisconsin cross-sections 
(Morrison and McCammon 1983) and relative elemental abundances as described by Anders and Ebihara (1982). For 
the spectra extracted from both regions we obtained acceptable fits using WABS 'x RAYMOND in both cases. In 
the case of the SeqID 50011000 observation, our fit values were kT = OAStoil keV and N u = 0.47±g;f| x 10 22 
cm -2 with x 2 /degrees of freedom = 42.6/149, while for the SeqID 57059000 observation our fit values were kT = 



0.571&1? keV and N R = 0.24lg;J| x 10 22 cm" 2 with x 2 /degrees of freedom = 71.6/149. We therefore conclude 
that the thermal emission from the northwestern rim of G156.2+5.7. may be adequately modeled using the a column 
density and a temperature corresponding to the average values for these two parameters derived from these two fits; 
that is, assuming iV H = 0.36 x 10 22 cm" 2 and kT = 0.52 keV. 

We next considered the RXTE PCA observation of the northwestern rim of G156.2+5.7: because of the large 
field of view of the PCA, this observation sampled flux from G156.2+5.7 as well as AX J0500+5238 and a sec- 
ond background galaxy, 3C 130. To properly account for emission from these two galaxies, we first extracted and 
modeled a GIS spectrum for AX J0500+5238 and fit the spectrum over the energy range of 0.7-6.0 keV using a 
power law and the WABS model for photoelectric absorption. We obtained an acceptable fit (x 2 /degrees of freedom 
= 90.5/448) to the spectrum with a photon index T = ISt 1 ^ and a column density N n = 0.57t°j7 x 10~ 22 cm" 2 
(both broadly consistent with Tomida 1999 and Yamauchi et al. 1999). In the case of 3C 130 (which was not sampled 
by the ASCA GIS observations), Hardcastle (1998) obtained an acceptable fit to the spectrum of this source using 
the RAYMOND model (kT = 2.9t\ keV) and a column density N H = 0.9+Q 2 x 10 22 cm" 2 , assuming a redshift of 
0.109. We adopt this model in our fitting procedure to describe the X-ray emission from 3C 130. 

To model the non-thermal component to the observed X-ray emission from the northwestern rim of G156.2+5.7, 
we considered two different models: a simple power law model and the SRCUT model (Reynolds 1998, Reynolds 
and Keohane 1999). This latter model describes a synchrotron spectrum from a power-law distribution (with an 
exponential cutoff) of electrons in a uniform magnetic field. The photon spectrum is itself a power-law, rolling off 
more slowly than exponential in photon energies. The SRCUT model can be used to estimate the maximum energy of 
electrons accelerated by the SNR: we assume that the relativistic electron energy spectrum N e (E) may be expressed 
as N c (E) = K E~ r e ~ E / EcutoH , where K is a normalization constant derived from the observed flux density of the 
region of the SNR at 1 GHz, V is defined as 2a+l (where a is the radio spectral index) and -E C utoff is tne maximum 
energy of the accelerated cosmic-ray electrons. A crucial advantage of this model is that a resulting fit may be 
compared with two observable properties of an SNR, namely its flux density at 1 GHz and a. Also, one of the fit 
parameters for this model is the cutoff frequency f cu toff °f tne electron synchrotron spectrum, defined as the frequency 
at which the flux has dropped by a factor of 10 from a straight power law. This frequency may be expressed as 



where B^q is the magnetic field strength of the SNR in fiG, assuming the electrons are moving perpendicular to the 
magnetic field. By using the value for f cu toff returned by SRCUT, an estimate for the maximum energy .E cu toff for 
the shock-accelerated electrons may be calculated. 

We therefore modeled the RXTE PCA spectrum using a four component model: a component for the thermal 
emission from the rim, a component for the emission from AX J0500+5238, a component for the emission from 3C 
130 and finally a component for the nonthermal emission from the rim. For the first three components, we adopt the 
previously-described models and parameter values, while for the final component we separately used and compared 
a power law model and an SRCUT model. Both of these models give acceptable fits to the non-thermal component 
of the spectrum: in the case of the power law model, for the photon index we obtain T = 2.0±g;| (x 2 /de grees of 
freedom = 16.2/15 = 1.08). Our value for the photon index is similar to values found for the bright rims of other 
SNRs that feature non-thermal spectra, such as G266.2-1.2 and G347.3-0.5 (Slane et al. 1999; 2001). This value 
is also comparable to the value obtained by Yamauchi et al. 1999 (T = 2.1) in their analysis of Ginga observations 
of G156.2+5.7. Using the SRCUT model also yields an acceptable fit (x 2 /degrees of freedom = 19.9/14 = 1.42) 
with values of 0.66±0.05 and 1.58±g ^ 

Jy for a and the normalization, respectively: both of these values are broadly 
consistent with the results of Reich et al. (1992). It is very remarkable that using the SRCUT model yields a fit which 
closely approximates the known radio properties of this SNR (note that the value for the normalization is almost half 
the total flux density at 1 GHz of G156.2+5.7 and reflects the bilateral radio morphology of this SNR). For the 
cutoff frequency parameter of this fit, we obtain ^ cut off = 2.41^0 23 x 10 17 Hz: assuming that the mean magnetic 
field strength of G156.2+5.7 is 14 fiG (Reich et al. 1992) and using Equation (1), we estimate that the maximum 
energy of cosmic-ray electrons accelerated along the northwestern rim of G156.2+5.7 is ss 32 TeV, well short of the 
"knee" feature of the cosmic-ray spectrum. In Figure 2, we present fits to the two ASCA GIS2+GIS3 spectra and the 
RXTE PCA spectrum that we have extracted for this analysis, using the SRCUT model to fit the non-thermal X-ray 
emission seen in the RXTE PCA spectrum. Additional X-ray observations of the northwestern rim of G156.2+5.7 




Hz, 



(1) 



Table 2. Gross Properties of Galactic SNRs with Non-Thermal-Dominated X-ray Spectra 



SNR 



SN 1006 



G28.6-0.1 



G156.2+5.7 G266.2-1.2 G347.3-0.5 



Other Names 
Distance (kpc) 
Age (yr) 



AX J1843.8— 0352 RX 



J04591+5147 RX J0852.0— 4622 RX J1713.7— 3946 



Size (arcmin) 

5*1.4 GHz (Jy) 
a (S u oc v~ a ) 

r 



2.2 1 
9971 

30 7 
19 7 
0.6 7 
2.1 10 
911 

«0.1 10 
la 1 



7 2 
<2700 2 
13x9 7 
3(?) 7 



2.1 2 




BQiG) 
n (cm -3 ) 
Progenitor Type 



<0.2 2 



12-16 8 - 150(^1°) 9 

«0.2 3 0.05+ «0.05-0.07 9 

II 4 II 5 



References: 1 - Winkler et al. 2002, 2 - Bamba et al. 2001, 3 - Yamauchi et al. 1999, 4 - Slane et al. 2001, 5 - Slane 
et al. 1999, 6 - Ellison et al. 2001, 7 - Green 2001, 8 - Reich et al. 1992, 9 - Pannuti et al. 2003, 10 - Allen et al. 
2001, 11 - Dyer et al. 2001. t - Calculated by Pannuti et al. 2003 using expressions derived by Slane et al. 2001 and 
assuming a distance of 1.5 kpc, a filling factor for a sphere of 1/4 and an explosion kinetic energy of 10 51 ergs. 

with improved sensitivity are required to more rigorously constrain its spectral properties. 

COMPARISONS WITH OTHER GALACTIC SNRS 

In Table 2, we present a comparison of the gross properties of five Galactic SNRs (SN 1006, G28.6— 0.1, 
G266.2— 1.2 and G347.3— 0.5, as well as G156.2+5.7) with X-ray spectra which feature non-thermal emission. All 
of the SNRs in this class are X-ray luminous but radio faint and are expanding into regions of low ambient density 
(n < 0.2 cm -3 ): it appears that searches for other SNRs in this class should concentrate on sources expanding into a 
low density media. However, we point out that two other SNRs which are known to possess non-thermal components 
to their X-ray spectra, Cas A and RCW 86, are expanding into high density ambient media. Therefore, the density of 
the ambient media is not the only factor in determining the strength of the non-thermal component to X-ray emission 
from SNRs. No clear evidence exists that any of these five SNRs are currently accelerating cosmic-ray electrons to the 
"knee" of the cosmic-ray spectrum: in fact, estimates for E cuto s for each of these SNRs often fall short by an order 
of magnitude or more (Reynolds and Keohane 1999). Additional observations of Galactic SNRs using the present 
generation of X-ray observatories are required in order to analyze the phenomena of non-thermal X-ray emission and 
cosmic-ray acceleration by SNRs. 
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